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The  use  of  polarization  diversity  process  requires  the  derivation  | 

i 

of  a mathematical  system  model  to  allow  evaluating  the  effects  of  this  | 

process,  employed  by  a jammer  or  by  a threatened  radar,  upon  jamming  | 

effectiveness.  The  system  model  consists  of  two  orthogonal  linear  j 

antennas  controlled  by  a random  process  in  amplitude  and  phase.  The 

! 

expected  value  of  the  received  power  due  to  interference  is  related  to 

i 

the  autocorrelation  of  the  received  signal.  The  received  power  is  ob-  ' 

tained  in  the  frequency  domain  from  the  spectral  behavior  of  the  total 

i 

polarization  diversity  processes  employed  by  both  the  jammer  and  radar  | 

and  the  spectral  characteristics  of  the  transmitting  and  receiving  ' 

systems.  The  spectral  behavior  of  tne  total  polarization  diversity 
process  is  constructed  by  the  convolution  between  the  polarization 
diversity  processes  of  both  the  jammer  and  radar.  The  spectrum  of  the 
polarization  diversity  is  spread  in  a convolution  maruier  by  the  trans- 
mitter spectrum.  This  convolved  spectrum  contributes  to  the  expected 
value  of  the  received  pov^er  only  within  the  bandwidth  of  the  receiving 
system.  It  is  concluded  that  for  several  typical  scenarios  polarization 
diversity  is  an  effective  jamming  technique. 
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THE  EFFECT  OF 


POLARIZATION  DIVERSITY 
ON  RECEIVING  SYSTEMS 
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I . Introduction 

Polarization  diversity  is  the  phenomenon  of  changing  the  polariza- 
tion characteristics  of  an  electromagnetic  wave.  The  results  of  this 
phenomenon  have  been  evaluated  and  implemented  primarily  in  radio  astron- 
omy with  some  applications  applied  to  radar  systems.  However,  the 
evaluation  of  the  effects  of  polarization  diversity  upon  reception,  pre- 
sented in  radio  astronomy  literature,  has  been  limited  to  the  case  where 
only  the  transmitted  polarization  is  varying  while  the  polarization 
state  of  the  receiving  antenna  is  kept  fixed.  The  objective  of  this 
thesis  is  to  derive  a mathematical  model  which  is  used  as  a tool  for 
evaluating  the  effects  of  polarization  diversity,  employed  by  a jammer 
or  by  a threatened  radar  system.  Once  the  mathematical  model  has  been 
derived,  it  is  shown  how  one  may  go  about  optimizing  the  likelihood  of 
reception  by  the  radar  system. 

Two  different  categories  of  polarization  diversity  are  presented. 

The  first  one  is  formed  implicitly  by  one  or  more  radar  systems  with 
different  fixed  polarization  states,  by  which  the  equivalent  polariza- 
tion state  seen  by  a jammer  becomes  a random  variable.  The  second  one 
is  when  a random  process  is  employed  explicitly  by  the  radar  system  in 
order  to  reduce  the  effectiveness  of  polarization  jamming.  For  both 
categories,  one  should  consider  the  statistical  behavior  of  the  polari- 
zation processes  when  evaluating  the  effects  of  those  processes  upon  the 
likelihood  of  reception. 
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The  derivation  of  the  time  domain  model  in  Chapter  II  is  done 
using  the  representation  of  a polarization  state  given  by  the  Stokes 
parameters  or  the  Poincare  sphere  and  using  the  theory  of  receiving 
partially  polarized  electromagnetic  waves  as  presented  the  radio 
astronomy  literature.  The  definition  of  the  Stokes  parameters  given 
in  the  radio  astronomy  literature  is  revised  to  enable  its  usage  in 
a more  general  situation  where  the  receiver,  as  well  as  the  transmitter, 
employs  a time  varied  polarization  diversity.  Implementation  consider- 
ations are  used  in  determining  the  system  model  used  for  polarization 
diversity.  The  general  transfer  function  of  the  receiving-transmitting 
system  model  in  the  time  domain  is  derived.  This  transfer  function 
is  reduced  to  the  function  given  in  the  radio  astronomy  literature 
once  some  assumptions  are  made.  Those  assumptions  are  listed  in  this 
chapter. 

In  Chapter  III,  the  system  model  derived  in  Chapter  II  is  trans- 
formed from  the  time  domain  representation  into  a statistical  represen- 
tation. The  averaging  process  which  is  used  for  evaluating  the  likeli- 
hood of  reception  is  done  by  the  ensemble  average  operation  rather  than 
by  a time  average  operation.  Since  the  expected  value  of  the  received 
power  is  related  to  the  autocorrelation  of  the  received  voltage,  the 
whole  system  is  described  in  a statistical  second  moment  representation. 
When  statistically  stationary  is  assumed,  the  Fourier  transform  of  this 
representation  leads  to  frequency  domain  representation.  The  effects 
of  the  polarization  diversity  processes  as  well  as  of  the  frequency 
characteristics  of  the  jammer  and  the  radar  systems  appear  in  the 
equation  of  the  average  received  power  in  the  frequency  domain.  The 
power  equation  is  simplified  for  various  classes  of  assumptions. 
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In  Chapter  IV,  several  feasible  scenarios  are  evaluated  with  respect 
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to  the  frequency  domain  equation  derived  in  Chapter  III.  Results  are 
obtained  for  several  cases  both  analytically  and  numerically. 
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II.  The  Time  Dependence  of  the 


Received  Power  on  Polarization  Diversity 

The  dependence  of  the  power,  received  by  a receiving  system,  on 
time  variation  polarization  diversity  which  is  employed  by  a transmitter, 
has  been  derived  in  the  radio  astronomy  literature  (Ref  i;  Ref  2).  This 
dependence  is  first  reviewed  in  this  chapter,  where  it  has  been  assumed 
that  the  receiving  system  has  no  variation  in  polarization.  The  depen- 
dence of  the  received  power  on  the  varying  polarization,  which  is 
employed  by  both  systems,  receiver  and  transmitter,  is  then  included. 

The  polarization  diversity,  discussed  in  the  radio  astronomy  liter- 
ature was  created  by  nature  and  there  has  been  no  need  to  discuss  the 
ways  to  implement  such  a phenomenon.  However,  since  this  study  is 
dealing  with  a man  made  polarization  diversity  system,  a feasible  model 
will  be  presented.  Thus,  the  derivation  of  the  effect  of  the  polariza- 
tion variations  will  be  related  to  that  model. 

Polarization  State 

A monochromatic  plane  wave  can  be  expressed  analytically  by  its 
electric  field  as 


E(r,t) 


where  Ic  is  the  propagation  vector  defined  by  k = kk. 


(1) 
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since  the  wave  is  plane  therefore  is  perpendicular  to  the  propaga-  ; 

tion  direction  k,  that  is,  = 0.  Without  loss  of  generality  k 

can  be  chosen  to  be  parallel  to  the  r axis  in  the  spherical  coordinate 
system  r,  0,  and  (ji . can  be  written  then  as 

E^  = §A  + $A2 

t 

= e|Ai|eJ“^  + $|A2|eJ“^  (2)  j 

I 

/V  /s 

where  0 and  <{)  are  the  unit  vectors  along  the  6 and  ({)  axes, 
respectively,  and  where  the  magnitudes  1ai|  and  IA2I,  as  well  as  the 
phases  Oj  and  0t2 , are  real  constants.  The  reason  for  choosing  the 
spherical  coordinate  system  is  explained  in  Appendix  A. 

It  follows  from  Equations  (1)  and  (2)  that  the  spherical  components 
of  E(r,t)  are  given  by  the  real  expressions 

Eg  = I Aj I cos(6  + ttj ) 

E^  = I A2 1 cos(o  + Oj ) 

i 

where 

6 = (i)t  - kr  (5)  i 

Eliminating  6 from  these  expressions,  one  obtains  the  following 
equation  (Ref  3:24-25) 


(3) 

(4) 


1 

E 

2 

E 

^ E 

E 

.1 

0 0 

d) 

A 

T 

lA  1 
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r 'j 

cosa  = sin  a 


(6) 


where 


a = Qj  - Oj 


(7) 


When  taking  Eg  and  as  coordinate  axes,  it  is  seen  that  Equation  (6) 

represents  an  ellipse  whose  center  is  located  at  the  origin  E„  = E.  = 0. 

D (p 

This  ellipse  is  illustrated  in  Figure  1.  The  polarization  state  of  the 
plane  wave  is  defined  by  this  ellipse  and  by  the  direction  of  the  rota- 
tion of  E(r)  along  the  ellipse.  The  direction  of  the  rotation  depends 
on  a.  The  ellipse  is  set  by  the  following  three  independent  parameters 
I Aj I , I Aj I , and  a . 

The  ellipse  can  be  represented  by  a different  set  of  parameters  which 
consists  of  the  semimajor  and  sem.iminor  axes  of  the  ellipse,  denoted  by 
a and  b respectively,  and  the  orientation  angle  ^ between  the  major 
axis  of  the  ellipse  and  the  9 axis  of  the  spherical  coordinate  system. 
This  set  is  illustrated  in  Figure  1.  The  relations  between  the  two  sets 
are  given  by  the  following  equations  (Ref  3:26-27) 


2 u2  ,2 

a + b = Ai  + 


Aa 


(8) 


and 


tan  2^  = 


2 All  1 

A2 

Ai 

2 

A2 

2 

cosa,  0 < < IT 


(9) 
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Figure  1.  Elliptic  Representation  of  Polarization  State 


which  were  derived  by  equating  the  two  different  representations  of  the 
ellipse. 

The  polarization  state  of  a plane  wave  may  be  defined  also  by  a set 
of  parameters  which  are  known  as  Stokes  parameters  (Ref  4).  The  Stokes 
parameters  of  a monochromatic  plane  wave  are  defined  by  the  following 
four  quantities. 


So  = 

+ 

|Aj*) 

1 

z 

o 

(10) 

Si  = (|aJ* 

- 

IMM 

1 

Z 

0 

(11) 

Sj  = (2|a,| 

lAj 

cosa) 

1 

Z 

o 

(12) 

Sj  = (2|aJ1 

|A,I 

sina) 

1 

Z 

0 

(13) 

where  is  the  intrinsic  impedance  of  an  isotropic  lossless  medium. 
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since  these  parameters  are  related  by  the  identity 


S ^ = Si^+  S2^+  83^ 

only  three  of  the  four  are  independent.  is  equal  to  the  Poynting 

vector  of  the  wave  since,  when  using  rms  quantities,  the  poynting 
vector  is 


_ _ _* 
S = E X II 


o 


(15) 


and  the  magnitude  will  be 


S = i(lAi|^  + UzP) 
^o 

= S 

o 


(16) 


The  Stokes  parameters  set  is  usually  normalized  such  that 


(17) 


and 


8 


12  _ 


z s 
o a 


aj  = 2"^  |Bi|  IB2I  cos(-e) 
o a 

83  = iBil  IB2I  sin(-6) 

o a 


(23) 


(24) 


(25) 


where  S is  the  magnitude  of  the  poynting  vector  or  the  power  density 
of  the  wave  radiated  by  the  antenna.  The  reason  for  choosing  a minus 
sign  in  front  of  S is  the  direction  of  the  received  wave  which  is 
opposite  to  the  direction  of  the  transmitted  wave  by  the  antenna,  by 
which  the  Stokes  parameters  have  been  defined. 

Alternatively,  the  Stokes  parameters  can  be  written  in  terms  of  the 
orientation  angle  'll  and  the  ellipticity  angle  X as  follows; 


Sj  = s^cos2Xcos2'iJ 


(26) 


s,  = s cos2Xsin2'l' 
^ o 


(27) 


S3  = s^sin2X  (28) 

Si , S2 , and  S3  can  be  interpreted  as  the  cartesian  coordinates  of  a 
point  on  a sphere  of  radius  s^  where  the  longitude  and  latitude  of  the 
point  are  2\p  and  2X  respectively.  The  ellipticity  angle  X is 
defined  by 

tanX  = ± 7 , -V4  < X < V4  (29) 

Si 
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where 


sin2X  = 


2IA1IIA2I 
I Ai  P+lA,  1 


1-2-  sina  = 


= s. 


(30) 


This  sphere  is  called  the  Poincare  sphere  and  it  is  illustrated  in 
Figure  2.  A point  on  the  sphere  represents  a polarization  state.  Linear 
polarization  corresponds  to  the  equator  of  the  sphere  or  to  zero 
ellipticity  angle  X*  Circular  polarization  corresponds  to  the  poles  of 
the  sphere  or  to  2X  = - '^/2  where  the  plus  sign  denotes  a right-handed 
circular  polarization  and  the  minus  sign  a left-handed  circular  polariza- 
tion. It  is  convenient  to  describe  the  polarization  state  by  either  one 
of  the  two  sets  of  angles  4*  and  X or  a and  Y where  Y is  defined 
as 


Y 


= arc tan 


0 < Y < T^/2 


The  relations  betvjeen  the  two  sets  are  given  by  (Ref  3:27) 


(31) 


tan2  4'  = tan2Ycosa 
sin2X=  sin2Ysina 


(32) 


The  two  sets  are  illustrated  in  Figure  2. 

In  summary,  the  polarization  state  has  been  defined  by  its  elliptic 
representation  and  by  two  parametric  representations,  Stokes  parameters 
and  Poincare  sphere.  The  representation  by  Stokes  parameters  is  used 
next  to  evaluate  the  received  power. 
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The  Effect  of  Polarization  Diversity  on  the  Reception  of  an  Antenna 

It  has  been  shown  by  Ko  (1962)  that  the  power  P available  from  an 
antenna  whose  effective  aperture  is  A^[a^],  when  a wave  of  polarization 
S[s^]  is  incident  upon  it,  is  given  by  (Ref  5) 

3 

P = y2SA^[^.][sJ  = y^SA^  I a.s.  (33) 

ex  1 ® i=0  ^ 

where  [a^]  is  the  transpose  of  [a^].  The  complete  derivation  of 
Equation  (33)  is  given  in  Appendix  A.  The  last  three  Stokes  parameters 
are  representing  the  three  components  of  a radius  vector  of  the  Poincare 
sphere.  Thus,  the  radius  vector  of  the  incident  wave  will  be  defined  as 

s = s^Sj  + SjSj  + SjSj  (34) 

12 


I 


and  the  radius  vector  of  the  antenna  will  be  defined  as 


/\  I 

3.  = + ^2^2  ^33^3  (35)  I 

1 


The  dot  product  of  the  two  vectors  could  be  replaced  by  the  cosine  of  the 


; angle  e between  the  vectors.  Thus 

I S'a  = Sjaj  + S2a2  + S3a3  = cose  (36) 

; 

i 

and  then 

P = J^SA^d  + cose)  (37) 

As  illustrated  in  Figure  3,  the  maximum  available  power  is  reached  when 
e = (})  or  2'n',  and  zero  power  is  received  when  e = it. 

A quasi-monochromatic  plane  wave,  which  consists  of  a superposition 
of  a large  number  of  statistically  independent  waves  with  a variety  of 
polarization  states,  is  said  to  be  partially  polarized.  The  frequency 
bandwidth  Af  of  such  a wave  is  assumed  to  be  very  small  compared  to 
the  mean  frequency  f.  The  electric  field  of  the  quasi-monochromatic 
plane  wave  may  be  uniquely  represented  by  the  analytic  form 

E(r,t)  = Re  E^(t)eJ^  -k-r) 

where  w denotes  the  average  value  of  the  frequency.  Since  the  band- 
width of  the  field  is  very  narrow,  will  be  a slowly  varying 

function  of  time  compared  to  Eg^t)  can  be  written  in  the  form 
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Figure  3.  The  Received  Power  versus  e 
E^(t)  = e|A2(t)|eJ“i^^^  + J|A2(t)|eJ“2^^^ 

where  the  magnitudes  )Aj(t)),  |A2(t))  and  the  phases  aj(t),  a 

are  slowly  varying  functions  of  time.  The  spherical  components  of 
are  given  by 

Eg  = |Aj(t)|cos(6  + ai(t)) 

E^  = |A2(t)|cos(6  + ai(t)  - a(t)) 

E = 0 
r 


where 


a(t)  = aj(t)  - a2(t) 
and  6(t)  = a)t  - k>r 


(39) 

(t) 

E(r,t) 

(40) 

(41) 

(42) 


(43) 


Although  the  magnitudes  and  phases  of  the  components  of  the  quasi- 
monochromatic  plane  wave  are  irregularly  varying  functions  of  time, 
certain  correlations  may  exist  among  them.  These  correlations  determine 
the  Stoke  parameters  and  consequently  the  polarization  state  of  the  wave. 
The  Stokes  parameters  of  the  quasi-monochromatic  plane  wave  have  been 
defined  as  the  time-averaged  quantities  (Ref  1) 


^ [<|Ai(t)l^>  + <lA2(t)l^>] 
^o 

Sj  = ^ [<lAj(t)|'>-<|A2(t)|'>] 
o 

Sj  = ^ <|Ai(t)| |A2(t)|cosa(t)> 

S3  = ^ <1 Aj (t) I I A^ (t) I sina(t)> 
o 


where 


<x(t)> 


T 

lim  1 
T^oo  2T 

-T 


x(t)dt 


(44) 

(45) 

(46) 

(47) 

(48) 


It  has  been  shown  that  the  Stokes  parameters  of  this  plane  wave  satisfy 
the  relation  (Ref  2:120) 

1 Si^+  S2^+  S3^  (49) 

By  dividing  through  by  S^,  the  normalized  Stokes  parameters  vector  will 
be  then 
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(50) 


F ’ 


i 


I 


[s^J  , i = 0, 1.2,3 

where  = 1.  The  degree  of  polarization  of  the  quasi-monochromatic 
plane  wave  is  defined  as 


s 

o 


= /sj^  S2^+  Sj^ 

(51) 

follows  that  for  a partially  polarized  wave 

be  represented  by 

the  Stokes  parameters  set 

s = 1 
o 

(32) 

Sj  = d cos2Xcos24' 

(53) 

Sj  = d cos2Xsin2'^ 

(54) 

s j = d sin2X 

(55) 

The  definition  of  the  Stokes  parameters  which  has  been  given  in  the 
radio  astronomy  literature,  can  be  applied  only  when  assuming  no  time 
variation  in  the  receiving  antenna  system  or  in  the  receiving  system 
itself.  However,  for  a complete  view,  one  should  consider  a general  case 
where  there  are  variations  in  time  in  the  receiving  antenna  system  and  in 
the  receiving  system  itself.  This  consideration  leads  to  the  definition 
of  the  Stokes  parameters  as  instantaneous  functions  of  time  while  the 
averaging  process  is  performed  on  the  instantaneous  power  function. 

Thus,  without  taking  into  consideration  the  processing  of  the  receiving 
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system,  the  average  received  power  will  be 


P = y,  <SA^[a.]ts.]> 

= '4  <SA  I a.s.>  (56) 

®i=0  ^ " 

where  the  averaging  process  will  be  defined  either  by  Equation  (48)  or  as 
a statistical  average,  depending  whether  the  variations  ai’e  deterministic 
or  random,  respectively.  The  averaging  process  ^^^ill  be  discussed  in 
detail  in  Chapter  III. 

When  defining  the  Stokes  parameters  by  the  9 and  <t>  components, 
there  is  no  consideration  how  one  would  go  about  implementing  the  plane 
wave  and  the  receiving  antenna  system  components  such  that  the  amplitudes, 
A^'s  and  B^'s,  can  be  controlled  independently.  The  necessity  of 
having  independent  control  of  the  amplitudes  results  from  the  objective 
of  this  study,  namely,  developing  an  optimized  polarization  div'ersity 
procedure.  In  the  following,  a feasible  implementation  of  a polarization 
diversity  model  will  be  considered. 


Polarization  Diversity  System  Model 

A wave  of  arbitrary  polarization  may  be  produced  by  combining  the 

I 

waves  radiated  by  a pair  of  crossed-dipole  antennas.  If  the  dipoles  are 
aligned  parallel  to  the  u and  v axes  of  an  arbitrary  coordinate 
system  (u,v,w),  then  the  components  of  the  wave  radiated  in  the  w 
direction  will  have  the  form 


[ 


Figure  4.  Polarization  Diversity  System  Model  ^ 

; 

in  which  |Aj | and  are  positive  real  magnitude  factors,  and  Oj 

and  ot2  are  the  phase  angles  associated  with  the  field  components.  The  i 

relative  magnitudes  and  the  phases  of  the  field  components  may  be  adjusted 

to  give  any  desired  polarization. 

The  system  which  will  be  used  for  evaluating  the  effects  of  polariza- 
tion diversity  upon  the  received  power  consists  of  two  orthogonal  linear 
antenna  channels  connected  to  a common  source  through  an  impedance 
matched  T adapter  which  is  called  a power  splitter.  This  system  will  be 
the  model  for  both  functions,  receiving  and  transmitting.  The  system  is  ! 

illustrated  in  Figure  4.  Each  antenna  channel  consists  of  a linear 

antenna  and  a magnitude  and  phase  controller  A.  which  is  controlling  ■ 

i 

the  input  or  output  current  of  the  antenna,  depending  whether  the  sys- 


The  characteristic  matrix  of  the  power  splitter  will  be 


C 0//T 

a/i/T 

= 

0//2  C 

0 

vT 

a//T  0 

C 

(59) 


where  r denotes  the  reflected  voltages  and  i denotes  the  incident 
voltages  at  the  associated  ports.  C denotes  the  reflection  coefficient 
and  O denotes  the  loss  coefficient.  When  the  splitter  is  matched  then 
^ = 0 and  0=1* 

The  effective  height  of  the  antenna  system  is  obtained  from 
Appendix  A as 


h = lAjje'j^ihj  + |Aj|eJ“2h2 


(60) 


(61) 


The  radiated  power  transmitted  by  this  antenna  system  is  derived  in 
Appendix  A and  is  given  by 

|hxf|*dj)  (62) 

n 

where  dJl  = sin0d6d(()  and  the  integration  is  done  over  a sphere. 
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Figure  5.  Equivalent  Circuit  of  a Transmitting  Antenna 

Due  to  matching  considerations  one  should  keep  the  total  radiating  resis- 
tance constant.  However,  since  the  radiating  resistance  of  this  model 
is  a function  of  the  A^'s,  there  has  to  be  a loss  resistance  which  will 
compensate  the  changes  in  the  radiating  resistance.  This  configuration 
is  illustrated  in  Figure  5.  Thus,  when  the  antenna  is  matched, 


R.  = R , 
1 rad 


(63) 


and 


(64) 


I 


Therefore 


"Md  ' ^<|A.I^|AJ') 


where  is  given  by  Equation  (A-37)  as 


■5^  Ihjxrl^df^ 


and  ^j^Qgg  represents  the  power  dissipated  in  the  antenna  system.  It 
has  been  assumed  that  the  t^vo  antennas  have  identical  input  impedances. 
One  can  verify  that  the  sum  of  the  radiated  powers  of  the  two  channels 
given  below  is  equal  to  the  radiated  power  of  the  total  system  assuming 
that  the  two  antennas  have  orthogonally  polarized  antenna  patterns. 


+ P4. 
ti  tz 


|I  I 2 

^o  2^^ 


\ f 

^ ^o  -*4^  Ih^x^l^dll 
, 


lb 

-Jx  OAjI^+IAjI^)  |hjXrj^dS7 
) , 


since 


I hxr  I 

n 


|(Ajhj  + A2h2)  X rl^df^ 


(AjhjXr  + A2h2xr) • (AjhjXr 


Ajhjxr)  df^ 


= A. 


IhjXrj^dQ  + Ia^I 


livrrdn 


(Iaj^  + 1aJ2)  ih,x;i^dJi 
n 


and 


(69) 


hjXr*h2xr  = 0 


(70) 


due  to  the  pattern  orthogonality. 

Assuming  the  same  model  for  both  the  receiving  and  the  transmit- 
ting systems,  one  obtains  the  total  system  as  illustrated  in  Figure  6. 

The  received  signal  v^(t)  will  be  formulated  by  using  Equation  (A-41) 
and  convolving  v^(t)  with  g^(t)  to  obtain 

Vo(t)  = g^(t)*[E^(t  - R/c)*h^(t)]  (71) 

where  h (t)  is  defined  as 
r 

h (t)  = B,(t)h  + B2(t)h  (72) 

r * ri  rz 


j. 
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E^(t)  is  obtained  from  Equations  (A-10)  and  (A-23)  by  considering 
in  the  frequency  domain  as 


E^(aJ)  = 
o 


_ + %)  -j(“  + %)R/c 

ATIH  t 


where  N^(a))  is  the  radiation  vector  in  frequency  domain  and  N^(ca)  = 

-rxrxN(ti))  while  N(a))  is  defined  in  Appendix  A by  Equation  (A-7). 

Since  u)<<a)  , 
o 


_ jiliiJi  e-0<"  ♦ N'CW) 

4TrR  ‘t^  ^ 


E^(t)  = 


- 5(t  _ B/c)*.N;(t) 

47IR  t 


e I.(t  - R/c)h^(t  - R/c) 


where 


Ii(t)  = v^(t)*g^(t) 


Due  to  the  pattern  orthogonality  one  can  approximate  h^(t)  by  h^  even 
for  angles  far  off  beam  center  for  many  such  antenna  types.  Thus, 


h^(t)  = Ai(t)h^^  + A2(t)h^, 


to  give 
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VQ(t)  = g^(t)* 

- ^ [v^(t-R/c)*g^(t)] 

h^(t-R/c) •h^(t) 

. 

(78) 


The  h .'s  and  the  h. .'s  are  assumed  to  be  constant,  with  respect  to 
ri  ti  > r 

frequency  and  time,  in  the  close  neighborhood  of  o)^.  By  using 
Equations  (72)  and  (77),  the  quantity  given  by 

h^(t-R/c) ■h^(t)  = [Aj (t-R/c)h^^  +A2 (t-R/c)h^^-  [B 1 ( t )h^j +B2 ( t )h^2 ^ 

= Aj (t-R/c)Bi (t) (h  'h  ) 
tin 

+ Aj (t-R/c)B2 (t) (h^j 'h^^ ) 

+ A, (t-R/c)Bi (t) (h  -H  ) 

^ * t2  ri 

+ Aj  (t-R/c)B2(t)(h^2’^r2^ 


The  received  power  is  given  in  Appendix  A by  Equation  (A-40)  as 


p (t)  - 
r^^'  " 4R, 


(80) 


As  illustrated  in  Figure  7,  the  maximum  received  power  will  be  achieved 
when 


R = R , + R- 
L rad  loss 


= R. 
1 
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(81) 


I 

4, 


Figure  7.  Equivalent  Circuit  of  a Receiving  Antenna 


where 


i 


Rrad  = (82: 

^loss=  ^i^^  ~ 


and  v(t)  is  measured  at  the  output  of  the  receiving  antenna  system. 

is  defined  by  Equation  (67).  Hence,  the  received  power  is  related 

* ' 
to  V (t)v  (t)  where 
o o 


V (t)v  (t)  = 
o o 


a)U 


4TrR 


|g^(t)*[  (v^(t-R/c)*g^(t))K^(t-R/c)  •E^(t)  ] (84 


This  is  the  general  transfer  function  of  the  received  power  of  the 
transmitting-receiving  system. 

It  can  be  shown  that  Equation  (84)  reduces  to  the  representation  of 
the  received  power  in  terms  of  the  Stokes  parameter  sets  of  the  trans- 
mitter and  the  receiver  as  a special  case. 
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In  the  radio  astronomy  literature  v^(t) 
g^(t)  was  assumed  to  be  C6(t)  where  C 
units  between  v^(t)  and  I^(t),  gp(t) 
antenna  was  assumed  to  be  lossless.  Thus 


was  assumed  to  be  constant, 
compensates  for  the  change  of 
was  assumed  to  be  5 (t),  and  the 


V (t)v  (t) 
= --4R-— 


rad 


1 

’ 

ojy 

rad 

47rR 

|h^(t-R/c) -h^Ct) I ^ 


In  the  case  where 


(85) 


\x= 


<j)h^^(e,(D) 


and  by  including  the  functions  of  0 and  (j)  in  the  respective  A^'s 

and  B . ' s 
1 

I h^(t-R/c) •h^(t) 1^  = I (t-R/c)Bj (t)  + Aj ( t-R/c iBj ( t ) | ^ 

= |a, (t-R/c)Bj (t)|2  + |A2(t-R/c)B2(t)|2 

+ 2Re[Aj (t-R/c)A*(t-R/c)Bj (t)B*(t) ] (86) 
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\ 


|h^(t-R/c).h^(t)r  = J^[|  Ai  (t-R/c)  |2+|A2(t-R/c)  |2][  |Bi(t)  |2+|B2(t)  |2] 

+ ;^[|Aj(t-R/c)|2-|A2(t-R/c)|2][|Bj(t)|2-|B2(t)H 

+ 2| Ai(t-R/c) I |A2(t-R/c) I |Bj(t) I |B2{t) Icos(a-B) 

(87) 

where 

cos(oi-S)  = cosacosB  + sinasing  (88) 


In  the  definition  of  the  Stokes  parameters  of  the  incident  wave,  B^ 
represents  the  electric  field  instead  of  the  effective  height 

Thus,  by  using  Equation  (A-23)  and  defining  the  field  as 


= AJ(tlh^_ 


(89) 


and 


the  received  power  will  be 


P 

r 


- — — Vi  ) S.a. 
4R  , 1 1 

rad  1=0 


3 

i 


(91) 


By  using  Equations  (A-47)  and  (A-49)  one  can  show  that  the  received 
power  will  have  the  form 


I 

[ 


P 

r 


s.a. 


1 1 


which  is  equivalent  to  Equation  (33). 


(92) 


i 

! 

1 
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III.  The  Average  Received  Power  as  a 
Function  of  the  Polarization  Diversity  Random  Processes 

The  system  which  was  illustrated  in  Figure  5 can  be  represented  as 
a transfer  system  where  v^(t),  h^(t) , and  h^(t)  are  inputs  and  v^{t) 
is  the  output.  This  system  is  illustrated  in  Figure  8.  The  space  trans- 
fer function,  g^Ct),  represents  the  attenuation  and  the  phase  changes 
in  the  signal  due  to  the  distance  R.  Thus 

G (0))  - - -j(a)+U)o)R/c  (93) 

’ ~ 47rR 

/ 

where  is  the  carrier  frequency, 


(aj+a3o)/c  = k + k 


Figure  8.  Transfer  Function  Representation 

As  it  was  stated  in  Chapter  II,  the  approach  for  evaluating  the  received 
power  will  be  by  averaging.  In  the  most  general  case  h^(t)  and  h^(t) 
are  assumed  to  be  random  processes.  One  could  evaluate  the  time  average 
received  power  by  assuming  ergodicity  and  evaluating  the  expected  value 
of  the  power.  In  any  case,  from  a statistical  point  of  view  one  may 
look  for  the  expected  value  of  the  power  since  it  will  tell  more  about 
the  likelihood  of  reception  at  any  time  than  would  the  time  average 
operation. 

The  received  power  per  unit  impedance  is  equal  to 

P (t)  = v^(t)v*(t)  (97) 

r o o 

and  the  expected  value 
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(98) 


f 


Pp  = E{v^(t)v^(t)} 


= P (0) 

V 

o 


where  R (t)  is  the  autocorrelation  of  v (t).  This  is  valid  when 
V o 

o 

assuming  stationary  processes  throughout  the  whole  system.  One  could 
represent  the  mean  power  in  terms  of  the  power  spectral  density  function 
as 


P = — 
r 217 


S (cj)da) 
o 


;99) 


where  S (oj)  is  the  Fourier  transform  of  R (t). 

V V 

o o 

It  is  desired  to  represent  S (w)  in  terms  of  the  spectral  density 

o 

functions  of  the  different  inputs  to  the  system,  such  that  one  could 
evaluate  the  mean  power  when  the  spectral  density  functions  of  the  inputs 
are  given.  As  shown  in  Appendix  B,  Equation  (B-42),  it  is  easy  to 
evaluate  the  spectral  density  function  of  an  output  to  a transfer  system 
as  a function  of  the  input.  Thus,  for  the  system  model 


^v  " ^v 
o r 


(100) 


and 


S (w)  = S ((jj)  |g  (oj)  I * 
i t 


(101) 
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However,  there  is  no  immediate  way  to  evaluate  the  spectral  density 
function  of  a multiplication  of  two  processes,  especially  when  they  are 


vectors.  The  spectral  density  function  of  v^(t)  will  be  derived  as 
follows.  Due  to  Equation  (A-41) 

= h^(t)-E^(t) 

= h^(t) • [E^(t)*g^(t) ] 

= h^{t) • [ (v^(t)h^(t) )*g^(t) ] (102) 

According  to  Equation  (96) 

V (t)  = h (t) "h  (t-R/c)v. (t-R/c)Q(R)  (103) 

r r t 1 

The  effective  heights  have  been  modeled  in  Chapter  II  as 

h, (t)  = A,(t)h,  + A,(t)h,  (104) 

t ^ t 1 ^ L 2 

and 


1 


A2(t)  = 

lA2(t) 

(107) 

B,(t)  = 

|B,(t)|e-J^i‘^^ 

(108) 

B2(t)  = 

lB2(t)  |e"J^2(t) 

(109) 

The  minus  signs  in  front  of  the  g^'s  are  due  to  the  opposite  directions 
of  the  receiving  and  transmitting  paths.  Thus 

V (t)  = Q(R)v. (t-R/c) [Aj(t-R/c)h  +A2(t-R/c)h  ] • [B j(t)h  +B^t)h  ] 

(110) 

and 

V (t-T)  = Q (R)v. (t-T-R/c) [A, (t-T-R/c)h,  +A2 (t-X-R/c )h,  ] 

r 1 1 ti  t2 

•[B*(t-T)h*^+B2(t-T)h^J  (111) 

The  resultant  autocorrelation  is 

R (T)  = e{v  (t)v*(t-T)} 

V r r 

r 

= E{h  (t) -h, (t-R/c)v. (t-R/c)Q(R) 
r t 1 

h^ ( t-T ) • h^  ( t-T-R/c ) ( t-T-R/c ) Q ( R ) } 

=|q(R)|^  E{h^(t)'h^  (t-R/c)h^  (t-T)'h^  (t-T-R/c) 

v^( t-R/c )v^  (t-T-R/c)}  (112) 
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It  is  reasonable  to  assume  that  v^(t)  is  statistically  independent 

from  h^(t)  and  h (t).  Thus 
t r 

R (t)  = 1q(R)1^R  (t)  E{h^(t) -h^Ct-R/clh^  (t-T)*h^  (t-T-R/c)} 

'^r  i 

(113) 


where  h . and  h i=l,2,  are  assumed  to  be  constant  in  frequency 
ti  ri 

for  (jj<<a)  and  to  be  random  variables  only.  The  quantities  A.  and 
o 1 

B . i=1.2,  are  assumed  to  be  random  processes  in  the  most  general  case. 

1 

As  a result,  the  h's  and  the  A’s  or  B's  are  statistically  independent. 
Therefore 


R (t)  = |q(R)1^R  (t)  E{h^(t) •h^{t-R/c)h^  ( t-x) "h^i t-T-R/c ) } 

'^r  i 

= (Q(R)|2r^  (t)  E{fAj(t-R/c)h^j->-A2(t-R/c)h^^]-[Bj(t)h^j 

i 

+ B2 ( t ) h^^ ] [ A* ( t-x-R/c ) h^*+A* ( t-x-R/c ) ] • [ B* ( t-T) h^* 

+ B*(t-X)h  *]}  (114) 

i r2 


or 


R (x) 

V 

r 


iQ(B)rB.  i i i i 

i=l  j=l  k=l  1=1 


V. 

1 


E{A.(t-R/c)Byt)A*(t-x-R/c)B*(t-x)(h^rH^  J(h^l^-h^^)*} 

. |g(R)l=»^  (T)  I Rijid)  Edhjj-K  (115) 

i i,j,k,l=l 
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where 


R'^hr)  = E{A.(t-R/c)A*(t-T-R/c)B.(t)B*(t-T)  } 

ik  1 R J -L 


(116) 


It  is  convenient  to  define  a cross-covariance  term  such  as 


ciJ^{T)  = E{[A.(t-R/c)A*(t-T-R/c)  - R^.j^(T)][Bj(t)B*(t-T)  - R^^^(t)]} 

(117) 


where 


Rtik(T)  = E{A^(t-R/c)A*(t-T-R/c)  } 


(118) 


and 


R (t)  = E{B.(t)B*(t-T)} 
rjl  J 1 


(119) 


B5'  this  definition,  can  be  represented  as  a function  of  the 

correlation  functions  of  the  receiver  and  transmitter,  R^^j^(t)  and 
R^^,^(T),  respectively.  Thus 


RJi(T)  . 


R.  (T)R  IT) 
tik  rjl 


(120) 


and 


R (t) 

V 

r 


|Q(R)1=R  (T)  I 

1 i,j,k,l=l 


(121) 
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or 


R (T)  = |q(R)1 

V 

r 


V. 

1 


-E  -h 


where 


i.j .k,l=l 


and 


i f J » ^ » 1 — 1 


dhx) 

ik 


E{(hj.-hrj)(hj^ 


V'*' 


When  transforming  into  frequency  space, 


where 


S (u)  = '?  [R„  (t)] 

V V 

r r 

S (u)  = ‘J'  [Ry  (t)  ] 
^'i 


and 
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(122) 


(123) 


(124) 


(125) 


(126) 

(127) 


%.h  <'•  * \-h 

t r t r t r 


= - ,T)J  - (T), 

t r t r 


= S„(CD)  + S (CD) 
R c 


(128) 


Substituting  Equation  (101),  Equation  (125)  becomes 


S (w)  1q(R)|^[S  (to)  Ig.  (w)  |2]*S^  ^ (CO) 

Vt  t h^-h^ 


(129) 


and 


P = 


Iq(r) r 


r 2n 


|G^(aj)  12{[s^  (oj)  |G^(a))  |2]*s-  (co)  }da) 

t t r 


(130) 


One  can  break  down  the  correlation  function  of  the  dot  product  h^'h^  by 
using  Equation  (B-27).  Thus, 


R^ik(T)  = E{A. (t-R/c)Ak(t-T-R/c) } 


■ itiV  * 


(131) 


and 


R^j.j(t)  = E{B^(t)B^(t-T)  } 


= 0 .T1  , + C .,  ( t) 
'rj  'rl  rjl' 


(132) 
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I 


where 


n^.  = E{A.(t-R/c)} 

= E{A|^(t-T-R/c)  } 

. e{b.(t)) 

and 

n*i  = E{B*(t-T)} 


(133) 

(134) 

(135) 

(136) 


Therefore 


* 1rj"rl 


(137) 


where  C..,(t)  and  C .,(t) 
tik  rj 1 

receiver,  respectively. 

In  summary, 


are  the  cross-covariance  of  transmitter  and 


P 

r 


S ((jj)d<jJ 
V 

o 


(138) 
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where 


(U))  = |G^(a))|2|G  |2{s^  .j-  (a))*[|G  (co)Ps  (o))]} 

o t r t 


and 


R^^(T)  = E{A^(t-R/c)A*(t-T-R/c)B^(t)B*(t-T)} 


In  a more  convenient  way, 


St-  ,7-  (u))  = S ((i))  + S (oj) 


where 


Sr((o) 


— « 4r 


with 


[E{  A^(t-R/c)A|^(t-T-R/c) }] 


[E{B  (t)B  (t-T)}] 

v)  ^ 


(139) 

(140) 


(141) 


(142) 


(143) 


(144) 


(145) 


f 


(146) 


ciJ(T)  = E{[A.  (t-R/c)A*(t-T-R/c)  - R^.|^(T)][Bj(t)B*{t-T)  - 


(147) 


Thus,  the  averaged  received  power  of  a receiving  system,  when  polarization 
diversity  processes  are  employed  by  the  receiver  and  the  transmitter, 
depends  upon  the  spectral  characteristic  of  the  correlation  between  the 
two  processes.  This  spectrum  is  spread  in  a convolution  manner  by  the 
transmitter  frequency  spectrum.  However,  this  convolved  spectrum  will 
contribute  to  the  average  power  only  within  the  bandwidth  of  the  receiving 
system.  The  parts  of  the  spectrum  which  are  outside  of  the  receiver  band- 
width are  wasted. 


Assumptions 

Equation  (138)  can  be  simplified  by  assuming  the  following 
assumptions : 

(1)  The  two  processes,  the  transmitting  and  the  receiving,  are 
statistically  independent. 

(2)  The  two  polarization  components  are  statistically  independent  for 
each  antenna  system,  both  transmitting  and  receiving. 

(3)  Only  one  of  the  systems  is  a random  process  and  the  other  is  a 
random  variable. 
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(4)  A typical  transmitter  where  v^(t)  = constant  and  a typical 

receiver  where  is  an  ideal  low  pass  filter  around  the 

carrier  frequency. 

(5)  S (t*))  is  a narrow  band  noise. 

'^t 

Assumption  (1)  yields 

= 0 (148) 


and 
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S.  ((jd)*S  .,(0))  E{(h^.-h  .)(h,  -h  .)*}  (149) 

tik  rjl  ti  rj  tk  rl 


where 


Stik(‘^)  = [E{A^(t-R/c)Aj^(t-T-R/c)}] 


(150) 


and 


Srji(“)  = [E{Bj(T)B*(t-T)}] 


(151) 


Since  h^  and  h^  are  statistically  independent  the  expected  value  of 
the  dot  product  of  their  components  can  be  simplified.  According  to  the 
model  which  has  been  used,  the  components  of  each  system,  the  transmitter 
and  the  receiver,  are  orthogonal.  Thus,  when  using  the  directions  of  the 


components  of  one  of  the  systems,  the  transmitter  for  example,  one  obtains 


^^^^ti’'^rj^^’^tk‘'^rl^  ^ '^^'^ti^rj^^  ^tk^rl^  ^ 


tk 


= E{h  h*  } E{h  . h*,  ) 

ti  tk  rj^.  rl^j^ 


(152) 


where 


h.  . = h . 'ti 

ti  ti 


(153) 


h . = h . 'ti 


^'Iti 


(154) 


ti  is  the  direction  of  and  the  other  terms  are  obtained  similarly. 


Thus, 


t r 


'ti  tk 


Due  to  Assumption  (2),  Equations  (131)  and  (132)  can  be  simplified, 
as  follows 


"ti"tk  * 

|n(J=  . c^i(T) 

II 

•H 

* 

\i\k 

, iA 

(156) 


where  C^^(t)  is  the  autocovariance  of  A^(t),  and 
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(157) 


Assumption  (3)  resolves  the  convolution  form  of  Equation  (149)  into 
a function  of  the  spectral  density  function  of  the  system  which  employs  a 
random  process  while  the  other  system  does  not  contribute  to  the  spectral 
characteristics  of  the  received  power,  except  for  a constant.  Since  one 
of  the  systems  behaves  as  a random  variable,  there  is  no  statistical 
dependence  between  the  two  systems  and  the  results  of  Assumption  (1)  can 
be  used  in  this  case.  Assuming  the  receiving  system  to  be  a random 
variable  yields  the  spectral  behavior 


^h  -h  J 

t r i,j,k,l=l 


S^.,  (cjj)  E{B  .B,  }E{h^.h. , }E{h  . h , 
tik  j 1 ti  tk  rj^.  rl^,^ 


1 


(160) 


Assumption  (4)  allows  an  explicit  insight  into  the  spectral  charac- 
teristics of  the  received  power.  This  insight  indicates  the  capability 
of  an  optimization  technique  which  is  discussed  in  detail  in  Chapter  IV. 
Since  = constant. 


(t)  = E{v^(t)v^( t-T) } 

= constant  (161) 


and 


S (w)  = D5(a)) 


(162) 


then  S may  be  written  as 
i 


(o))  = D5(a))  iG^(aj)  1^ 
i 


(163) 


45 


with  the  receiver  definition 


lG^(a))|2 


0 , |u)l>w 


one  can  resolve  Equation  (130)  to  obtain 


(164) 


^ |Q(R)|"|G^IMciJO)12  ' (aj)da) 
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(165) 


^h  *h  I ^<S(a))  ] = 3j^  ((jj-co* ) 1g^((jo' ) [ ^6(co' )dai' 

"t  "^r  ■*  "t  r 


= %-h 

t r 


(166) 


The  results  of  Assumption  (4)  are  changed  only  slightly  with  Assump- 
tion (5).  The  only  difference  is  that  the  spectral  density  function 

^h  ‘h  spread,  in  a convolution  mariner,  by  the  spectmim  of  the 

t r 

noise  of  the  transmitter. 
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IV.  statement  of  Problems 


One  should  consider  two  opposite  problems  concerning  polarization 
diversity.  Both  problems  are  dealing  with  the  question  of  how  to 
optimize  the  likelihood  of  reception  of  a receiving  system  when  a 
polarization  diversity  phenomenon  is  employed  by  both  the  receiving  and 
the  transmitting  systems.  The  first  problem  is  how  to  maximize  the 
likelihood  of  reception  of  the  receiving  system  in  order  to  effectively 
disrupt  reception.  This  problem  is  referred  to  as  the  Electromagnetic 
Counter  Measure  (ECM)  problem.  The  second  problem  is  how  to  minimize 
the  likelihood  of  reception  of  interference  by  the  receiving  system 
in  order  to  avoid  jamming.  This  problem  is  referreo  to  as  the  Electro- 
magnetic Counter-Counter  Measure  (ECCM)  problem.  In  both  cases,  there 
are  two  opposing  systems,  the  receiver  and  the  jammer,  where  one  system 
is  assumed  to  have  some  kind  of  polarization  diversity  technique  and 
the  designer  has  to  find  an  optimized  polarization  diversity  technique 
for  the  opposing  system  such  that  the  likelihood  of  reception  is 
optimized  in  the  desired  sense.  The  "assumed"  polarization  diversity 
technique  can  employ  a random  process  in  the  most  general  case,  or  it 
can  be  simply  a random  variable,  meaning  that  the  polarization  can  take 
on  different  values  with  each  value  having  some  probability  of 
occurrence . 

Optimization  Factors 

One  of  the  objectives  of  this  study  is  to  point  out  some  optimization 
techniques  by  employing  random  processes  with  the  polarization  diversity 
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Figure  9.  Maximum  Spectrum 

system  model  which  vas  developed  in  Chapter  II.  One  can  use  the  result 
of  the  most  general  case  as  it  appears  in  Equation  (137).  However,  for 
the  sake  of  simplicity,  several  assumptions  are  made  in  order  to  point 
out  an  optimized  technique  of  polarization  diversity.  Some  specific 
limitations  on  implemetation  are  also  considered.  Using  Assumption  (4) 
to  define  a simple  system,  one  obtains  the  spectral  function  of  the 
received  power  as  given  by  Equation  (164) 


<ii 

c 

"r  = " f \.h 

^ t r 
-0) 
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f 

(166) 


where  S—  (w)  is  defined  by  Equation  (139). 

^t  "r 

From  the  ECM  point  of  view,  optimal  reception  occurs  when  most  of 

the  energy'  of  the  spectral  density  function  -h  between  -<ii 

t r 

and  as  illustrated  in  Figure  9.  Optimal  reception,  from  the  ECCM 


J 
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Figure  10.  Minimum  Spectrum 

viewpoint,  is  achieved  by  spreading  the  spectrum  of  the  received  power 


such  that  most  of  its  energy  lies  outside  the  bandpass  to  cu^, 

as  illustrated  in  Figure  10.  The  behavior  of  the  amplitude  of  (oj) 

”t  "r 


is  another  factor  in  the  optimization  procedure.  To  minimize  the  likeli- 
hood of  reception,  one  may  desire  the  ajnplitude  of  (oj)  to 

"t  "r 

approach  zero.  This  is  achieved  when  the  two  processes, 

h^(t),  are  statistically  orthogonal.  Maximum  likelihood  of  reception 

occurs  when  the  two  processes  are  statistically  dependent.  Obviously, 

it  is  difficult  to  implement  a complete  dependence  and  one  may  assume, 

for  simplicity,  that  the  two  processes, 

statistically  independent  as  stated  by  Assumption  (1). 


Evaluating  Tvpical  Scenarios 


Some  typical  scenarios  are  evaluated  according  to  the  system  model 


which  was  presented  in  the  last  two  chapters.  Four  types  of  scenarios 
are  considered. 


(a)  The  pol '^rization  state  of  the  jammer  is  fixed  while  the 


polarization  state  of  the  receiver  is  assumed  to  be  a random 
variable . 

(b)  The  configuration  is  the  same  as  in  the  first  case  except 
that  it  is  considered  as  an  ECCM  problem. 

(c)  The  jEimmer  employs  a random  process  in  order  to  overcome  the 
ECCM  technique  used  in  the  second  case. 

(d)  The  receiver  employs  also  a random  process  i.n  order  to  reduce 
the  jamming  effectiveness. 

For  evaluating  the  received  power  in  the  first  scenario  one  may  go 
back  to  the  representation  of  the  polarization  state  by  the  Poincare 
sphere,  or  specifically  to  Equation  (37),  where  the  power  is  given  as 
a function  of  the  angle  between  the  polarization  states,  of  the  jammer 
and  the  radar,  such  that 

P = EiVid  + cose)}  (167) 

where  P is  the  normalized  average  power.  Five  different  cases  of 
scenario  (a)  are  listed  in  Table  I.  As  seen  in  the  table,  in  most  of 
the  cases,  the  likelihood  of  reception  is  optimized  from  the  jammer 
point  of  view.  It  may  be  expected  that  the  polarization  state  of  the 
radar  will  be  adjusted  to  be  orthogonal  to  the  jammer  such  that  the 
average  power  is  reduced.  This  status  is  categorized  under  scenario  (b). 

The  obvious  action  that  may  be  taken  for  improving  jamming  effec- 
tiveness is  employing  a random  process.  This  case,  listed  as 


Circular  Circular  1 , 
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right  hand  is  equal  to  0.5 
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Figure  11.  Random  Binary  Transmission 

scenario  (c),  is  evaluated  analytically  and  numerically  where  both 
results  are  presented  and  compared.  The  scenario  is  defined  by  the  four 
amplitudes  of  the  system  model  such  that 


Aj  (t) 

A2(t) 

Bj(t) 


B2(t) 


1 

1 

0 


(168) 


where 


ai(t)  = {-Tr/2,  + ir/2} 


(169) 


and  it  is  a pseudo  random  process  called  random  binary  transmission  as 
illustrated  in  Figure  11.  The  B's  are  selected  as  such  for  the  ease  of 

evaluation,  h is  assumed  to  be  a random  variable  with  equal 

rl 
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Figure  12.  Autocorrelation  of  Random  Binary  Transmission 


Probability  of  being  oriented  in  G or  ({)  . h and  h are  assumed 

Ll  L2 

A ^ 

to  be  oriented  in  9 and  i , respectively.  Thus, 


Rjl<T,  . Xtl  . R^^(T) 
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where 


R (t)  = 


1 - , Kl  < T 


(171 


T > T 


as  illustrated  in  Figure  12.  The  spectral  density  function  of  R (t), 


S . (oj) , is  given  as 


c*  / /.•  \ _ 'r 


sin*  (iTTf ) 


Figure  13.  Power  Spectral  Density  of  Aj(t) 


as  illustrated  in  Figure  13.  Thus,  the  average  power,  P^,  is  given  a. 


P = '4 
r 


= Vz 


[6(f)  + (f)]df 
A 

1 


(173) 


assuming  the  spectrum  of  the  jammer  to  be  within  the  radar  bandwidth. 

Most  of  the  energy  of  the  spectrum  S (f)  is  in  the  range  lf|5 

Aj  T 

Hence,  when  the  radar  is  band  limited,  jamming  is  effective  only  when 


T > 


(174) 


where  f^  is  the  cutoff  frequency  of  the  radar  bandwidth. 

Scenario  (c)  is  also  evaluated  numerically  by  the  use  of  SAP3 
program  found  in  the  AFIT  computer  library.  The  random  process  employed 
by  the  jammer  is  constructed  to  be  similar  to  the  Random  Binary  Trans- 
mission illustrated  in  Figure  11.  The  process  is  based  on  the  sign  of 


54 


a zero  mean,  unit  variance,  Gaussian  random  noise  filtered  by  several 
types  of  low  pass  filters.  Figure  14  illustrated  the  process,  created 
by  a 320KHZ  filter.  Figures  15  and  16  illustrate  the  time  autocorrela- 
tion and  the  spectral  density  functions  of  this  process.  Assuming 
erogdicity,  time  and  ensemble  autocorrelation  may  be  equated  which 
allows  comparison  of  the  analytical  and  numerical  analysis.  The  spec- 
tral density  function  is  integrated  along  the  frequency  domain,  as 
illustrated  in  Figure  17.  The  value  of  the  integral  for  a certain 
frequency  is  equivalent  to  the  expected  received  power  obtained  by 
passing  the  process  through  the  receiver  with  a bandwidth  equal  to  that 
frequency.  Figure  18  illustrates  the  spectral  density  function  con- 
structed by  filtering  by  640KHz.  Figure  19  illustrates  the  expected 
received  power  versus  the  receiver  bandwidth.  In  Figures  20  and  21  the 
spectral  density  function  and  the  received  power  are  illustrated  as 
before  except  for  a process  created  by  a 1.28MHz  filter.  The  filtering 
of  the  process  has  the  same  influence  on  jamming  effectiveness  as  the 
changing  of  T in  the  analytical  analysis  does.  Thus,  for  a certain 
receiver  bandwidth  the  expected  received  power  decreases  when  the  band- 
width of  the  random  process  increases.  Thus,  as  the  bandwidth  of  the 
polarization  diversity  process  becomes  larger  than  the  bandwidth  of  the 
radar  system,  the  jamming  effectiveness  goes  down. 

The  last  scenario,  (d),  is  analyzed  numerically  due  to  its  complex- 
ity. The  jam.ming  process  is  assumed  to  be  as  before  and  the  radar  polar- 
ization diversity  process  is  assumed  to  be  similar  to  the  jamming 
process.  The  radar  process  is  assumed  to  be  created  by  a 640KHz  filter 
for  all  cases.  The  process  of  the  radar  system  is  employed  by  the 
polarization  ai.iplitude  controller  Bj(t).  In  Figures  22  through  27,  the 
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joint  spectral  density  functions  and  the  expected  value  of  the  received 
power  are  given  respectively  for  all  three  cases  of  jammer  filtering 
values,  320KHZ,  640KHz,  and  1.28MHz.  For  each  case,  the  spectrum  of 
the  jammer  process  is  convolved  with  the  spectrum  of  the  radar  process 
to  obtain  the  joint  spectrum.  The  joint  spectrums  are  spread  due  to 
the  convolution.  Figures  22,  24,  and  26  demonstrate  the  gradual  spread 
of  the  joint  spectrum  due  to  the  increase  of  the  spectral  bandwidth  of 
the  jammer  process.  Figures  23,  25,  and  27  illustrate  the  increasing 
dependence  of  the  expected  value  of  the  received  power  upon  the  receiver 
bandwidth  due  to  the  increase  of  the  spectral  bandwidth  of  the  jammer 
process.  Thus,  as  was  concluded  in  scenario  (c),  jamming  effectiveness 
is  reduced  as  the  bandwidth  of  the  jammer  process  becomes  larger  than 
the  bandwidth  of  the  receiver.  However,  the  radar  process  as  well  as 
the  jammer  process  dominate  jamming  effectiveness  since  both  are 
playing  the  same  role  in  creating  the  spectral  behavior  of  the  joint 
polarization  diversity  process. 


V.  Conclusion 


The  use  of  a polarization  diversity  process  requires  the  deriva- 
tion of  a mathematical  tool  which  will  allow  one  to  evaluate  the 
effects  of  this  process  upon  the  received  power.  Partial  derivation 
exists  in  the  radio  astronomy  literature  for  a special  case  where  only 
the  transmitter  employs  a random  process  to  produce  the  polarization 
diversity.  Since  random  processes  are  assumed  to  be  used  for  polari- 
zation diversity,  the  expected  value  of  the  received  power  has  to  be 
a function  of  the  statistical  characteristics  of  the  processes.  The 
time  average  operation  used  in  radio  astronomy  is  not  as  indicative 
of  the  statistical  characteristics  of  the  processes  at  an  instant  of 
time  as  much  as  may  be  the  ensemble  average. 

Any  polarization  state  may  be  produced  by  combining  the  waves 
radiated  by  a pair  of  cross-dipole  antennas.  This  knowledge  was 
used  to  construct  a system  model  of  polarization  diversity.  The 
mathematical  dependence  of  the  received  power  upon  polarization  diver- 
sity was  derived  by  the  use  of  this  system  model.  The  general  trans- 
fer function  of  the  received  power  was  introduced  first  in  the  time 
domain.  This  function  included  the  behaviors  of  the  transmitter  and 
the  receiver  by  using  their  time  domain  transfer  functions.  The  trans- 
fer function  of  the  received  power  was  transformed  into  the  frequency 
domain  where  the  received  power  was  measured  by  its  expected  value. 

The  likelihood  of  reception  was  introduced  as  a function  of  the  spec- 
tral behavior  of  the  polarization  diversity  process  and  the  spectral 
characteristics  of  the  transmitter  and  the  receiver.  The  spectral 
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characteristic  of  the  total  polarization  diversity  process  is  construc- 
ted by  the  convolution  between  the  polarization  diversity  processes  of 
both  the  transmitter  and  receiver  systems.  The  frequency  spectrum  of 
the  polarization  diversity  is  spread  in  a convolution  manner  by  the 
transmitter  fi-equency  spectrum.  The  convolved  frequency  spectrum  con- 
tributes to  the  expected  value  of  the  received  power  only  within  the 
bandwidth  of  the  receiving  system. 

Two  opposite  approaches  exist  in  optimizing  the  likelihood  of 
reception  of  interference,  maximization  and  minimization  depending 
whether  the  ECM  or  the  ECCM  problem  is  considered,  respectively.  The 
optimizations  factors  of  the  polarization  diversity  process  are  first 
the  distribution  of  spectral  behavior  of  the  process  with  respect  to 
the  receiver  bandwidth  and  second  the  amplitude  of  the  spectral 
behavior. 

Four  feasible  types  of  scenarios  were  evaluated  with  respect  to 
the  system  model  of  polarization  diversity.  The  first  two  scenarios, 
consisting  of  random  variables,  were  actually  resolved  by  the  use  of 
the  Poincare  sphere  representation  of  polarization  states.  Results  of 
the  last  two  scenarios,  involving  random  processes,  were  obtained 
analytically  and  numerically.  In  both  methods  of  analysis  it  was  shown 
that  by  controlling  the  rate  of  change  or  equivalently  the  spectral 
behavior  of  the  random  processes  employed  by  both  the  transmitter  and 
receiver  polarization  diversity  systems  one  could  achieve  either  maxi- 
mization or  minimization  of  the  likelihood  of  reception. 

There  is  a definite  disadvantage  to  the  radar  system  if  a random 
variable  type  of  polarization  diversity  technique  is  employed  by  a 
jammer  when  the  polarization  state  of  the  radar  system  is  fixed.  If 
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the  polariation  state  of  the  radar  system  is  controlled  to  nullify 
jammer  interference  by  polarization  orthogonality,  the  polarization 
diversity  technique  employed  by  the  jammer  is  even  more  crucial  and 
should  be  chosen  to  be  a random  process  rather  than  a random  variable. 
This  essentially  guarantees  optimal  likelihood  of  reception  by  the 
radar.  In  the  case  where  the  polarization  state  of  the  radar  is  a 
random  process  used  to  reduce  jamming  effectiveness  in  the  manner  of 
spread  spectrum,  the  action  that  should  be  taken  by  the  jammer  is  to 
correlate  its  process  to  the  radar  process.  This  action  compresses 
the  spectral  behavior  of  the  total  polarization  diversity  process  in 
a manner  similar  to  the  spectral  compression  of  the  radar  signal. 

This  thesis  did  not  take  into  consideration  the  problems  one 
might  face  when  choosing  a polarization  diversity  process  for  either 
the  jamming  or  the  receiving  systems.  These  problems  are  to  retain  the 
effectiveness  of  existing  ECM  or  ECf'M  techniques  while  optimizing  the 
likelihood  of  reception  of  the  jamming  signal  by  the  radar  and  to 
develop  a simple  and  feasible  system  to  actually  accomplish  the  polar- 
ization diversity.  These  problems  are  suggested  for  further  study. 
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Appendix  A 


Derivation  of  Equation  33 


The  purpose  of  this  appendix  is  to  show  the  detailed  derivation  of 
Equation  33.  Some  brief  derivations  have  been  written  in  the  literature 
(Ref  1).  Considering  the  vector  potential  at  a point  p,  with  a distance 
R from  the  origin  of  a radiating  system,  as  illustrated  in  Figure  A-1, 
one  obtains  the  equation 


T .-jkr" 


(A-1) 


where  A is  the  phasor  of  the  vector  potential  A such  that 


A = AeJ"’'^ 


(A-2) 


In  the  far  field,  when  R»r',  the  following  assumptions  can  be  made.  For 
amplitude  consideration  the  assumption  is 


r"  = R 


(A-3) 


but  for  phase  consideration 


r"  = R-r'cos 


(A-4) 
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Figure  A-1.  Geometry  of  a Radiating  System 


where 


cosip  = cos6cos6'  + sinBsin©'  cos((('-<J)' ) 


(A-S) 


Thus 


A 


e r -T-  ^jkr'cos'l' 

4TTR  J ^a^ 


(A-6) 


The  integral  is  a function  of  the  antenna  configuration,  current  distribu 
tion,  and  direction  of  the  current.  This  integral  is  defined  as  the 
radiation  vector  N such  that 


N = 


-r  jkr'cosi/' 


( A — 7 ^ 


V 


Thus 


A = ® 


-jkR 


471 R 


(A-8) 


It  is  convenient  to  describe  the  radiation  vector  as  a multiplication 
of  the  input  current  and  an  effective  height  vector  of  the  antenna 

h.  For  a thin  linear  antenna,  this  becomes 


-I 


— jkr ' cosi|J  , , ,,  , 
1 e'^  ds'dl' 

a 


1'  s' 


i (s')ds' 
s 


^(l,)^jkr'cos.]^dl' 


{A-9) 


or 


N = I. 
1 


y(l,)^jkr'COSl|7 


= I.h 
1 


(A-10) 


where  s'  is  a surface  intersecting  the  feed,  f(l')  is  the  current 
behavior  along  the  length. 


h = 


^(l,)^jkr'cos4. 


(A-11) 


and 


A = " 


-jkR 


I.h 


47rR  i 


(A-12) 
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In  a source-free  medium,  the  plane  wave  form  of  the  far  field  is 


it  . H .-j’'"'' 
o 


E = E 

o 


(A-13) 


(A-14) 


where  k = kr  and  H and  E are  constant  vectors.  Thus, 

o o 


E = — Vx  H 
Jtoe 


(A-15) 


and  when  neglecting  terms  which  decrease  faster  than  1/R  in  the  far  field 


E ==  — ^-jkxH) 


(A-16) 


Since 


H = V X A 


{A-17) 


and 


A = A e 
o 


therefore 


-jk*r 


(A-18) 


H = -jkxA 


(A-19) 


and 
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E = - — r X r X A 

juje 


(A-20) 


. ^ ^ II  , 

E = jwy  r X r X A 


= -jwy  A' 


(A-211 


where 


A'  - AqO  + A^v*) 


(A-22) 


Replacing  A'  by  h'  according  to  Equation  (A-12)  leads  to 


4^R  1 


(A-23) 


where 


h'  =-rxrxh 


he  6 + h^<f, 


with  no  r component.  Thus 


(A-24) 


E = Ege  + E^<{, 


(A-25) 


where 


. ^ I.h 

4trR  1 0 


(A-26) 


I 

i . . 


and 


(A-27) 


(A-28) 


E _ _ iilH  I h 

~ 4TTR  14 


The  poynting  vector  of  the  electric  field  has  the  form 


p = Re[E  X H*] 


assuming  rms  magnitudes.  The  magnitude  of  the  vector  poynting  is 


p = ^ [E-E*] 
o 


since  the  magnetic  field  is  perpendicular  to  the  electric  field  in  a plane 
wave.  Thus  the  magnitude  of  the  poynting  vector  of  an  electromagnetic 
field  produced  by  an  antenna  will  be 


■>-f  S'ii 

o '• 


where  h is  the  effective  height  of  the  antenna.  Since 


f-  = 0)41? 

o 


(A-30) 


(A-31) 


p becomes 


80 


2 


P = z. 


■2AR 


(h* -h'*) 


= Z 


{hi 

Hr 


(A-32) 


Hence,  the  total  power  radiated  by  an  antenna  will  be 


P = 


p*ds 


IT  277 


=1 


pR^dfJ 


o o 


where 


(A-33) 


dfi  = sin6d6d(}) 


(A-34) 


Thus 


P = Z 


fll  11 

2 ^ 

r 2 

f f 

' i' 

2A 

1 

0 0 

(h'-h'*)df2 


(A-35) 


The  power  radiated  by  the  antenna  can  be  expressed  also  by  its  feed 
current  and  effective  radiation  resistance.  Thus,  for  rms  magnitudes 


rad 


{A-36) 


assuming  no  losses  in  the  antenna,  therefore  the  effective  radiation 
resistance  will  be 


I 
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Figure  A-2.  Equivalent  Circuit  of  an  Antenna 


rad 


IT  2ir 

I {h’-h'*)dj^ 
o o 


(A-37 


The  maximum  power  received  by  an  antenna  is  reached  when  its  internal 
impedance  is  equal  to  the  conjugate  impedance  of  the  radiation  impedance 
or  that 


R , = R. 
rad  1 


and 


X j = 
rad 


(A-31 


The  model  circuit  of  an  antenna  is  illustrated  in  Figure  A-2.  Thus,  the 
maximum  power  received  by  the  antenna  will  be 


i 


i 

I 

r 


p 

r 


* 

V V 
PC  PC 

4R  j 
rad 


(A-40) 


When  the  antenna  is  net  matched  the  ppwer  will  differ  within  a cpnstant 
pnly.  Sinclair  has  shewn  by  using  the  reciprecity  theerem  that  the  epen 
circuit  vcltage  of  a receiving  antenna  can  be  expressed  as  (Ref  6) 


V = E ’h 
PC  e 


(A-41) 


where  h is  the  effective  height  vecter  cf  the  receiving  antenna  and 
is  the  incident  electric  field  upen  the  antenna.  Thus,  substituting  v^^ 
in  Equaticn  (A-40)  yields 


P 

r 


(E  •h)(E  -h)* 
p p 


4R 


rad 


(A-42) 


Since  in  far  field  there  is  nc  r cempenent  cf  the  incident  field 


P 

r 


4Rrad  ^^9^® 


(A-43) 


The  directive  gain  cf  an  antenna  is  expressed  by 


D(9,(J)) 


_ _* 

E*E 


J-  f 

417 


(E'E*)df2 


(A-44) 


If  the  antenna  is  100%  efficient,  which  means  nc  mismatch  cr  ccnductive 
Icsses,  the  gain  and  the  directive  gain  are  the  same.  If  the  antenna  is 
net  100%  efficient  then  the  gain  will  be  Icwer  than  the  directive  gain  by 
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I 


a factor  which  depends  on  the  efficiency.  Assuming  a perfect  antenna, 
the  expression  of  the  gain  becomes  then 


_ _* 

E*E 


(E*r*')dfi 

Q 


(A-45) 


Since  h'  is  relative  to  E by  Equation  (A-23),  therefore  the  gain  can  be 
represented  as  a function  of  h' 


0(0,4)) 


h' -h'* 


^1: 

(h'-h'*)dn 

(A-46) 


Substituting  Equation  (A-46)  in  Equation  (A-37)  yields  a simpler  expression 

for  R ,,  that  is 
rad 


4Trz  ^ ^ I * 
o h' •h' 

rad  “ 4X®  0(0  ,(j)) 


(A-47) 


Substituting  Equation  (A-47)  in  Equation  (A-43)  yields 


P 

r 


x^G(e.4))  7 

4Tr  o 


(E  •h)(E  -h)* 
o o 

(h' -h'*) 


(A-48) 


By  multiplying  the  numerator  and  the  denominator  by  the  power  density  of 
the  plane  wave  at  the  receiving  antenna  as  in  Equation  (A-29)  one  obtains 
as  a result  the  received  power  as 


which  is  the  effective  aperture  of  the  antenna  and  S is  the  real  part  of 
the  poynting  vector  of  the  plane  wave  or  the  first  Stokes  parameter  as  in 
Equation  (16).  The  components  of  the  field  and  the  effective  height  can 
be  expressed  as  a magnitude  and  phase  functions  of  time 


(A-52) 


(A-53) 


(A-54) 


= B2(t) 


= |B2(t)  (A-55) 

where  the  minus  sign  in  front  of  the  B^(t)'s  is  due  to  the  opposite 
direction  of  the  transmitting  path  from  the  receiving  path.  Thus,  by 
substituting  Equation  (A-52)  through  (A-55)  in  Equation  (A-50)  one 
obtains 

A S 

" (|Aj(t)|2+|A2(t)|2)(|Bi(t)i2+|B2(t)i2)  |Ai(t)|"|Bi(t)|2 

+|A2(t)l^|B2(t)|2 
+ |Ai(t)  I |A2(t)  I |Bj(t)  I |B2(t) 

+ |Ai(t)  I |A^(t)  I |Bi(t)  I lB2(t) 

(A-56) 


or  by  reordering  the  terms, 


( |Aj(t)  l'=+|A2(t)  |2)(  |Bj(t)  |2+|B2(t)  |2)  l'ilAi(t)  |2  |Bj(t)  1^ 

+ )i|Aj(t)  |2|B2(t)  |2+J^|A2(t)  |2|Bj(t)  |2+>i|A2(t)  |2|B2(t)  |2 
+ )S|Aj(t)  |2|Bj(t)  12-y2|Aj(t)  |2lB^(t)  |2-Vj|A^(t)  |2|Bj(t)  [2 
+ y2\k^{t)\^\E^(t)  |2+)$4Re  [Ajt)  I |A2(t)  I |b  j(t)  I |B2(t)  ^ 

(A-57) 


where 


a(t)  = aj(t)  - a^(t)  (A-58) 

and 


6(t)  = B^(t)  - B2(t)  (A-59) 

Since 

= costa(t)  - 3{t)] 

= cos  a(t)  cos  6(t)  + sin  a(t)  sin  3(t)  (A-60) 


and  by  substituting  Equations  (10)  through  (13)  and  Equations  (22)  through 
(25)  into  Equation  (A-57)  one  obtains 

Sa  + S,a,  + S,a,+  Sa 

r,  V.  c O O 11  2 2 3 3 


= y^A  S[a.  ][s.  ] 
e 1 1 


i = 0, ... ,3 


(A-61) 
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Appendix  B 


Definitions  and  Equations  in  Random  Fields 

This  appendix  provides  definitions  and  equations  in  random  fields 
for  the  statistical  terms  used  throughout  the  thesis.  When  further  infor- 
mation considering  these  terms  is  needed  one  may  refer  to  Papoulis  (Ref  7). 

Random  Variables 

A variable  which  is  a function  of  the  result  of  a statistical  exper- 
iment, in  which  each  outcome  has  a definite  probability  of  occurrence,  is 
called  a random  variable. 

Probability  Distribution  Function.  The  probability  distribution 
function,  F (x),  of  a random  variable  X is  defined  as  the  probability 

A 

that  X < X,  that  is 


F^{x)  = P(X<x) 


(B-1) 


In  case  of  an  experiment  with  two  random  variables,  X and  Y,  the  joint 
probability  distribution  function  is  defined  as 


= P(XOc,  Y<y) 


(B-2) 


Probability  Density  Function. 


random  variable  X,  f (x)  or  f(x) 

\ 


fj(x)  = 


dx 


The  probability  density  function  of  a 
for  simplicity,  is  defined  as 


(B-3) 
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The  joint  probability  density  function  is  defined  as 


1 1 


I 

1 


t- 


= 3x3y 

Statistical  Dependency . Two  random  variables,  X and  Y,  are  statis- 
tically independent  if,  for  any  x and  y 

P(X<x,  Y<y)  = P(X<x)  P(Y<y)  (B-5) 


or 


(B-6) 


If  the  two  random  variables  are  statistically  dependent,  the  joint  proba- 
bility density  function  will  be 


(B-7) 


where  ^xjy  conditional  probability  density  function  of  X given 

that  Y is  known  and  is  defined  as 


X.  i..\  lin>  P[x-Ax<Xs:x  given  Y=y] 
= Ax^ 


and  similarly  for  f^iY^yl^)* 

y ; X 


(B-8) 


i. 


i 

i 


89 


t ; 


Expected  Value.  The  expected  value  E{x},  or  the  mean  n.  of  a 
random  variable  X,  is  defined  as 


E{X}  = 


X f^(x)dx 


(B-9) 


If  X can  have  only  discrete  values  x then 

n 


e{x}  = y X p 
^ n'^n 
n 


(B-10) 


where  p^  is  the  respective  probability  of  the  value  x^.  The  expected 
value  of  a function  of  a random  variable,  Y=g(X),  is  given  by 


E{Y} 


* 

g(x)f  (x)dx 

J X 


(B-11) 


or 


E{Y}  = I g(x^)p^ 


(B-12) 


The  expected  value  of  a function  of  two  random  variables,  g(X,Y),  will  be 


E{8(X,Y)}  = 


I g(x,y)f^Y(^»y)dxdy 


(B-13) 


where  f^(x,y)  is  defined  by  Equations  (B-6)  or  (B-7).  The  expected 
value  of  a linear  combination  of  N random  variables  is  equal  to  the 
same  linear  combination  of  their  expected  values.  Thus, 
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I 


r ■ 


N N 

{ I a.X.}=  I a.E{X.} 
i=l  i=l 


(B-14) 


Uncorrelated  Variables.  If  X and  Y are  uncorrelated  then 


e{xy}  = e{x}  e{y} 


(B-15) 


However,  when  g(X,Y)  = u(X)v(Y)  and  X and  Y are  uncorrelated,  it 
does  not  necessarily  follow  that  u and  v are  uncorrelated.  But  if 
X and  Y are  statistically  independent  then 


E{g(X,Y)}  = E{u{X)}  E{v(Y)} 


(B-16) 


Orthogonality . Two  random  variables,  X and  Y,  are  called  orthogonal 


if 


E{XY}  = 0 


(B-17) 


Variance.  The  variance  a*  of  the  random  variable  X,  is  defined  by 


= E{(x-n)*} 


(x-n)  f(x)dx 


(B-18) 


or 
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I 


0^  = I (x^-n)* 

n 

when  X has  only  discrete  values  x^.  The  variance  can  be  expressed  by 

0^  = e{x}  - e2{x} 

= e{x^}  - (B-20) 

Random  Processes 

A random  process  is  defined  as  a random  variable  in  which  to  each 
outcome  of  the  experiment  a time  function  is  assigned. 

Stationary  Processes.  A random  process  is  stationary  in  a strict 
sense  if  its  statistics  are  not  affected  by  a shift  in  the  time  origin. 
This  means  that  the  two  processes  X(t)  and  X(t+G)  have  the  same 
statistics  for  any  e.  The  two  processes,  X(t)  and  y(t),  are  jointly 
stationary  if  the  joint  statistics  of  X(t),  Y(t)  are  the  same  as  the 
joint  statistics  of  X(t+e),  Y(t+e)  for  any  e.  The  expected  value  of  a 
stationary  process  X(t)  is  a constant, 

E{x(t)}  = n 

= constant  (B-21) 

Correlation.  For  the  purpose  of  this  thesis,  the  autocorrelation  of 
a stationary,  complex  process  X(t)  is  defined  by 

Rjj(T)  = E{X(t+T)  X*(t)}  (B-22) 


1 

1 

(B-19)  ! 
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assuming  the  real  and  the  imaginary  parts  of  X(t)  to  have  the  same 
statistics.  The  cross-correlation  of  two  stationary  processes,  X(t) 
and  Y(t),  is  defined  as 


= E£x(t+T)  Y*(t)} 


(B-23) 


A process  X(t)  is  stationary  in  a wide  sense  if  its  expected  value  is 
a constant  and  its  autocorrelation  depends  only  on  x = ti  - tj . Two 
processes  are  jointly  stationary  in  a wide  sense  if  each  one  of  them 
satisfies  Equations  (B-21)  and  (B-22)  and  their  cross-correlation  depends 
only  on  T = tj  - t^,  as  in  Equation  (B-23).  In  the  following  discussion, 
only  stationary  processes  will  be  considered.  Two  processes  X(t)  and 
Y(t)  are  called  uncorrelated  if,  for  any  x 


Covariance.  The  autocovariance  of  a process  X(t)  is  defined  as 


C(T)  = E{(x(t+T)  -n)(X*(t)  - n*)}  (B-25) 

The  cross-covariance  of  two  processes  X(t)  and  Y(t)  is  defined  as 


C^(t)  = E{(X(tH-x)  -nv)(Y*(t)  - n*)}  (B-26) 


I 


p 


If  X(t)  and  Y(t)  are  uncorrelated  then,  for  any  x, 


(B-28) 


Cjj(T)  . 0 


(B-29) 


If  X(t)  and  Y(t)  are  orthogonal  then 


"xy<''>  = 0 


{B-30) 


Ergodicity.  X(t)  is  ergodic  in  the  most  general  form  if  its  time 
averages  are  equal  to  the  ensemble  averages  (i.e.,  expected  values). 

Time  Invariant  System.  Transformation  of  a process  X(t)  by  a time 
invariant  system  without  memory  leaves  the  statistics  of  the  results  Y(t) 
similar  to  those  of  X(t).  The  output  of  such  a system  will  be  stationary 
only  if  the  input  is  stationary  in  the  strict  sense.  The  expected  value 
of  Y(t)  = g[X(t)]  will  be 


w 

E{Y(t)}  = I g(x)  f^(x;t)dx 


and  the  autocorrelation  will  be 


Ry(T)  = E{Y(t+T)  Y (t)} 


»||  g(x, 


) g(x  ) f(x  ,x  ;T)dx  dx 
2 12  12 


(B-31) 


(B-32) 
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where  and  are  dummy  variables,  f(Xj,X2;T)  is  the  joint 

density  function  defined  as 


f(Xj,X2;T)  = 


3 F(Xj,X2;T) 


(Jx  dx. 
1 2 


(B-33) 


F(x^,x^;T)  = P{x(t+T)  < Xj,  X(t)  i-x^} 


(B-34) 


Power  Spectral  Density  Function 

The  power  spectral  density  S(a))  of  a process  X(t)  is  the  Fourier 
transform  of  its  autocorrelation,  R(T): 


= 1 


R(  T)dT 


(B-35) 


S(w)  is  a real  and  nonegative  function.  From  the  Fourier  inversion 
formula  follows  that  R(  T)  can  be  expressed  in  terms  of  S(u)  by 


= H f 


S(a))e'^'**^  du) 


(B-36) 


With  T=0,  the  above  yields 


VI 

■ 

2Tr 


S(a))dw  = R(0) 


= E{X(t)X  (t)} 


(B-37) 


This  is  equal  to  the  "average  power"  of  the  process  X(t).  The  cross-power 
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spectral  density  function  S (uj)  of  two  processes  X(t)  and  Y(t)  is 

A I 

the  Fourier  transform  of  their  cross-correlation: 


S„(u.)  . 


The  Fourier  inversion  formula  gives 


XY' 


= 1 


= I dw 


XY 


and  with  T=0, 


2TT 


= E{x(t)Y  (t)  } 


(B-38 


(B-39 


(B-40 


If  the  processes  X(t)  and  Y(t)  are  orthogonal,  then,  due  to  Equation 
(B-30) 


Sxy(‘*>)  = 0 (B-41 

Linear  Systems.  The  power  spectral  density  of  the  output  of 

a linear  system  with  system  function  H(jii))  is  given  by 

SyCu))  = S^(a))  |H(ja))P  (B-42 

where  S (oj)  is  the  power  spectral  density  of  the  input.  The  linear 
X 

system  is  illustrated  in  Figure  (B-1). 


Sx(a)) 


H(  ju)) 


-Y(t) 

= s (u))|H(ja))|* 


E{Z(t)Z*(t-T)}=  E{X(t)Y(t)X*(t-T)Y*(t-T)} 

= E{X(t)X*(t-T)Y(t)Y*(t-T)} 

= E{X(t)X*(t-T)}E{Y(t)Y*(t-T)}+C2,(T) 


= Rx(t)Ry(t)  + C^CT) 


(B-43) 


where  ^2^^^  called  the  cross  covariance  and  will  be  defined  as 


C^Cr)  = E{[X(t)X  (t-T)-Rx(T)][Y(t)Y  (t-T)-RY(T) ] } 


= E{X(t)X  (t-T)Y(t)Y*(t-T)}  - Rx(T)Ry(T) 


(B-44) 


The  product  is  illustrated  in  Figure  (B-2).  If  X(t)  and  Y(t)  are 


Figure  B-1.  The  Linear  System 


Multiplication  of  Two  Processes . The  autocorrelation  of  the  product 
Z(t)  of  two  random  processes,  X(t)  and  Y(t),  can  be  expressed  as 


statistically  independent  then 


C2(T)  = 0 


(B-45) 
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X(t)  ^Z(t)  = X(t)Y(t) 


Figure  B-2.  The  Product  of  Two  Random  Processes 


From  Equation  (B-27), 


= C^(T)  + 


and 


Ry(T)  = C^lT)  + iTlyl^ 

Thus 

E{X(t)Y(t)X*(t-T)Y*(t-T)}  = R^(T)Ry(T)  + C^{T) 

= [Cj^(T)  + |njjr][CY(T)+|nYl^]  + c^{T) 
= c^(t)Cy(t)  + InYl^c^(T) 


(B-46) 


(B-47) 


+ 


|n„rc„(T)  + + c,(t)  (b-48) 
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